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Abstract. A regularized a— system of the Nonlinear Schrodinger Equation (NLS) with 2<x nonlinear 
power in dimension N is studied. We prove existence and uniqueness of local solution in the case 
1 < a < N 4 _2 an d existence and uniqueness of global solution in the case 

1 < f < jf- When a -> 0+, 

this regularized system will converge to the classical NLS in the appropriate range. In particular, the 
purpose of this numerical rcgularization is to shed light on the profile of the blow up solutions of the 
original Nonlinear Schrodinger Equation in the range j? < a < A, and in particular for the critical 
case a = 
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1. Introduction 

The Nonlinear Schrodinger equation (NLS): 

iv t + Av + \v\ 2a v = 0, i-el", tel, (1) 

v(0) = v , 

where v is a complex-valued function in x R, arises in various physical contexts describing wave 
propagation in nonlinear media (see, e.g., [14], [21], [22] and [26]). For example, when a = 1, equation 
(1) describes propagation of a laser beam in a nonlinear optical medium whose index of refraction is 
proportional to the wave intensity. Also, the Nonlinear Schrodinger Equation successfully models other 
wave phenomena such as water waves at the free surface of an ideal fluid as well as plasma waves. In all 
cases, it is interesting to note that Eq. (1) describes wave propagation in nonhomogeneous linear media 
with self- induced potential given by \v\ 2a . 

As it is mentioned above, the a = 1 case is particularly interesting for laser beam propagation in 
optical Kerr media. Depending on the dimensionality of the space upon which the beam is propagating 
in, the wave dynamics can be either "simple" or "intricate" . In one space dimension, the NLS equation 
is known to be integrable and possesses soliton solution that preserves their structure upon collision [1]. 
The picture in two-dimensional (2D) space is totally different. The 2D NLS equation is not integrable, 
hence no exact soliton solutions are known. Instead, the 2D NLS equation admits the waveguide solution 
(also known as Townes soliton) v(x,y,t) = i?(r) exp(ii) with r = \J x 1 + y 2 where R > satisfies the 
nonlinear boundary value problem 

+ -^-R + R 3 = 0, ^(0)=0, lim R(r)=0. (2) 
ar r ar ar r^+oo 
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Importantly, the L 2 norm (or power in optics) of the Townes soliton defines a critical value for blow up. 
If initially the beam's power is larger than that of the Townes soliton, ||i>o||f,2 > H-RHI2, then the beam 
undergoes a finite time blow up. If on the other hand ||i>o||r2 < | |? 2 3 then the wave will diffract. 
Various mechanisms to arrest collapse have been suggested such as nonparaxiality [6], or higher order 
dispersion [7]. 

As a result, an important issue that arises in the mathematical study of the NLS is the question 
of local and global existence of solutions, their uniqueness, as well as the profile of blow up solutions. 
Knowing answers to such questions may have some consequences on possible physical observations of 
phenomenon governed by the NLS and in validating its derivation. In the works of Ginibre and Velo 
[10] and Weinstein [24], it is proved that equation (1) has a unique global solution when < a < j?, 
and that it has a unique global solution for "small" initial data for the critical case a = -w. The 
proof of global existence uses the fact that the energy Af(v) = J RN \v(x,t)\ 2 dx and the Hamiltonian 

H(v) = J RN ^|Vu(x,t)| 2 — '"^^ J dx are conserved quantities of the dynamics of (1). In the case of 

a > jji Glassey [11] proved that there exist solutions that develop singularities in finite time. In recent 
years there was an intensive remarkable computational work concerning the blow up for the critical case 
a = jj. For instance, Merle and Raphael have obtained a sharp lower bound on the blow up rate for 
the L 2 norm of the NLS in M. N (see [15] and references therein). Moreover, Fibich and Merle [8] studied 
self-focusing in bounded domains using a combination of rigorous, asymptotic and numerical results. 

Instead of the potential \v\ 2a , physicists consider self-gravational potential (see, e.g., [18], [20]) and 
come to a new system: Schrodinger- Newton equation (SN): 

iv t + Av + ipv = 0, lei", iel, (3) 

-a 2 At/; = \v\ 2 , 
v(0) = Vq, 

where a > is a real constant. System (3) is a Hamiltonian system with the corresponding Hamiltonian 
Tl{v) — J RN ^|Vu(x,t)| 2 — V^^H^M)! ^ x anc [ can ^ obtained formally by the variational princi- 
ple = where v* denotes the complex conjugate of v. System (3), or at least its stationary 
state, has been studied [17], [23], [19]. This coupled system of equations consists of the Schrodinger 
equation for a wave function v moving in a potential ip, where ip is obtained by solving the Poisson 
equation with source p = \v\ 2 . It can be thought of as the Schrodinger equation for a particle moving 
in its own gravatational field. As in the NLS, the energy J\f(v) = J RN \v(x, t)\ 2 dx and the Hamiltonian 

Ti.(v) = J RN ^|Vu(x,i)| 2 — iiSiOMsSl ^ d x are a i so conserved in this system. The question of existence 
and uniqueness of local and global solutions for system (3) has not been answered completely yet. 

Inspired by the a— models of turbulence (see, e.g., [3], [4], [5], [9], [12], [13] and references therein), 
we introduce a generalization of (3), the Schrodingcr-Hclmholtz (SH) regularization of the classical NLS: 

iv t + Av + u|u| a " 1 w = 0, xel" teM., (4) 

u-a 2 Au = M CT+1 , 

V(0) = Vq, 

where a > and a > 1. System (4) is a Hamiltonian system with the corresponding Hamiltonian 
(x,t)\ 2 - u ^ v a ^ +1 ) dx and can be obtained formally by the variational principle 

= ) where again v* denotes the complex conjugate of v. In this system, we can regard the 

wave function v moves in a potential u|u| cr_1 , where u is obtained by solving the Hchnholtz elliptic 
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problem u — a 2 Au = \v\ a+1 . Observe that the energy Af(v) = J RN \v(x,t)\ 2 dx and the Hamiltonian 

= Jr« (jVv(.T, i)| 2 — "( X '*)M>'*)I J ^ x are conserved in this system. When a = 1, we have the 

potential u as in the SN with the only difference that the Possion equation is modified as a Helmholtz 
equation. So we consider this system as a generalized system of SN. A more important fact is that when 
a = 0, one recovers the classical NLS, therefore we regard this system as a regularization of the classical 
NLS. In this paper we focus on the case a > and in our subsequential works, we will investigate the 
behavior when a — > + . In particular, we will investigate the case a = j^, which is not completely 
understood. 

In this paper, we will study the question of local and global existence of unique solution for system 
(4). Specifically, we will prove the short time existence of unique solution, when 1 < a < jfz^, ( we define 
once and for all N _ 2 = oo when N < 2). Moreover, we will show global existence of unique solution 
when 1 < a < The proof will follow the ideas of [10] and [24] and use the important fact of the 
conservation of the corresponding energy and the Hamiltonian of (4). All the proofs presented here will 
apply directly to system (3) as well. So simultaneously we have the same results for system (3): we have 
short time existence of unique solution when 1 < a < N 4 _ 2 , and global existence of unique solution when 
1 < a < jf- Comparing to the results of the classical NLS (1) (a < for local existence and a < 
for global existence), one expects these "better" results for (3) and (4) since the nonlinear terms in (3) 
and (4) are milder than that of the classical NLS (1). The parametre a plays an important role in our 
proofs. In a subsequential paper, we will investigate numerically the blow up profiles of the NLS, in the 
relevant range of er, when a — > + . 

In section 2, wc will introduce some essential notations and definitions, and some preliminary results 
that will be used throughout the paper. Following the work of Ginibre and Velo [10], we prove in section 
3 local (in time) existence and uniqueness of solution for system (4) using the contraction mapping prin- 
ciple. In section 4, we will extend the local solution to global existence, for 1 < a < after establishing 
the required a priori estimates for the H 1 norm of the solution, which remains finite for every finite 
interval of time. 

2. Notations and Preliminaries 

In this section we introduce some preliminary results and the basic notations and definitions that will 
be used throughout this paper. 

We denote by || • || p the norm in the space LP = L P (R N ) (1 < p < oo), except for p = 2 where the 
subscript 2 will be omitted. We will denote by {■, ■) the scalar product in L? . The conjugate pair p, p' 
satisfies the relation i + ^ = 1. For any real number I, wc denote by H l = H l (M. N ). the usual Sobolev 
space. Of special interest is the H 1 Sobolev space with the norm defined by 

11*= / (i + \e)\m\ 2 dt (5) 

or equivelently, 

IMI^HMP + Iiv^l 2 . (6) 

We denote by HitH^.j, = fe\ a \<k Jwn \D a u\ p dx) ^ P , 1 < p < oo, for u belongs to the Sobolev space 
W k,p (M. N ). For any interval / of the real line R, and for any Banach space B, we denote by C(I,B) 
(respectively C&(/, £>)) the space of continuous (repectively bounded continuous) functions from / into B. 

In this paper C and C a will denote constants which might depend on various parameters of the problem. 



4 



Y.CAO, Z. H. MUSSLIMANI, AND E.S. TITI 



They might vary in value from one time to another, but they are independent of the solution. When 
it is relevant we will comment on the asymptotic behavior of these constants as they depend on the 
corresponding parameters. 

First, we recall some classical Gagliardo-Nircnberg and Sobolcv inequalities (see, e.g., [2]). 
Proposition 1. (1) For any TV > 1, we have 

\\v\\ q <C\\v\\ 1 -^ N \\Vv\\ S ^ N for every v G H 1 , < ^— -TV < 1 (7) 

Lq 

\\v\\ q < C\\v\\ W 2, m for every v G W 2 ' m , q > m, 2m > TV (8) 

HI, < C\\v\\ w2 , m for every v e W 2 ' m , - >-- ^ > 0, q < oo (9) 

q m TV 

In particular, (10) 
h\\ q < C\\v\\ W 2,* = C\\v\\ H 2 for every v G H 2 , 2 < q < oo, TV < 3. (11) 

(2) For TV < 2, 

\\v\\ q < C\\v || ff i for every v G H 1 , 2 < q < oo. (12) 

With these inequalities at hand, we can process the nonlinear term. Let us rewrite the term 

f(v) = u^^v = B(\v\" +1 )\v\' r - 1 v, (13) 

where B = (I — a 2 A) _1 , the inverse of the Hclmholtz operator. Then / is a locally Lipschitz mapping 
from H 1 into L r , for some r G (2, -^35], where £ + X = 1. 

Proposition 2. Let TV > 1 and 1 < a < N *_ 2 . For every v±,v 2 G H 1 C If, where r depends on 

the given a and belongs to the range r G (2, ™ 2 \ (we consider 2N ^ 2 as 00 when TV < 2), we have 

\\f(vi) - f(v2)\\r> < k\\vi - v 2 \\ r , where k = C a {\\vx\\ H i + ||^2 1| jj-i ) 2cr and i + -K = 1, for some constant 
C a . 

Before we prove this proposition, we will state the following Lemmas: 

Lemma 3. Let TV > 1 and 1 < a < N A _ 2 . For every i>i,v 2 ,v G H 1 C L r , where r depends on the given 
a and belongs to the range r G (2, j^ ], we have 

\\B(\v\' y )\vr 1 v(v 1 -v 2 )\\ r , <C a \\v\\%4v 1 -v 2 \\r, (14) 
||B(|«r- 1 «)|«r- 1 «(«i-«a)||r' < CJvf^Wv, - v 2 \\ r . (15) 

Proof. First, denote 

h = \\B(\vn\vr i v(v 1 -v 2 )\\ r ,. 

Case 1. TV < 2: 

By Holder's inequality, we have 

h < \\B(\vn\vr i v\\ r , Pl \\v 1 -v 2 \\ rl , 1 

= \\B(\vn\vr 1 v\\^\\v 1 ~v 2 \\ r , (16) 
where in the last equality, we choose 71 = 1 — 1 > 1 and + — j- = 1 such that r'ji = r and r' (3\ = ^35- 

By Cauchy-Schwarz inequality, we have 

mwM^vw^ < \\B(\vn\\ - iimi - . (17) 
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Now, for the elliptic equation u — a 2 Au = / in Mr, we have the regularity property [16], [25] 

||«||w2,p < C a ||/||p for any 1 < p < oo, (18) 

where C a depends on N,p and a, and C a ~ -\ as a — ► + . Moreover, for a fixed, C a ~ p as p — » oo. 
Since ^ > 2 > by (11) and (18), we have 

iib(ioii_- < ciiB(it,nn^, 2 
< c a \\\vn 
= c a \\v\r 2a , 

and 

IIMl^HMI^- 

Since a > 1, combining the above two terms and applying (12), we have 

h <C a || \\ Vl -v 2 \\ r . 

Case 2. N > 3: 

Applying Holder's inequality, we have 

h < \\B(\vn\\r' B A\\vr 1 v\\ r ' A\vi-V2\\r' 71 

= ||B(H' r )||r'» 1 ||»||^ / 9 1 ||«l-^||r' 71 , 

where j- + j- + ± = 1. 

Now by (8), (9) and (18), we have 

IISODIU < C\\B(\vn\\ w2 , m 
< CyiMlm 

= I w 1 1 am ' 

where we require -4- > - — 4f when - — 2. > o, or r'#i > m when - — 4r < 0, for m > 1 to be 
determined later. 

Therefore, we obtain 

Ii<C a \\v\^ m \\v\\^ 1 \\v 1 -v 2 \\ rl71 . (19) 
Now, by requiring am = err' Pi = r ,r y\ = r, we have 

r-1 

> 1 a < r - 2 



r - a - 2 



r~ 1 

/3i = >l=>cr<r-l 

(T 

71 = r — 1 > 1 
r 

771 = — > 1 
(7 

AT + 2 1 r - 2 

o" < r — 1 or <7 > . 

2N 2 

Since 2 < r < j^z^, we conclude that a < jyz^, i.e., 

*1 <CaH|l?i||«X-«a||r. 

By exactly the same steps, inequality (15) follows readily. □ 
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Lemma 4. Let N > 1 and 1 < a < jrio ■ F° r ever U vi,v 2 ,v e H 1 C L r , where r depends on the given 
a and belongs to the range r G (2, we have 

WBdvr+^vr'im-i^Wr, <C a |H|^i||«i-«a|| r , (20) 
^(M^H^V^ - v 2 )\\ r , < C a \\v\\% v 2 \\ r . (21) 

Proof. Denote 

I 2 = \\B(\v\° +1 )\vr 1 (v 1 -v 2 )\\ r ,. 

Case 1. N < 2: 

By Holder's inequality, we obtain 

h < \\B(\v\^ 1 )\vr- 1 \\ r , 01 \\v l -V 2 \\ r ^ 

= \\B(\vr +i )\vr i \\^\\v 1 -v 2 \\ r , 

where we choose the same /3i,7i as in (16). 

Now, when a > 1, again by Holder's inequality, we have 

IIBCK+^r 1 !!^ < \\B{\v\^)\\^ 2 \\\v\°-^ 12 

(22) 

= ||S(|«r +1 )||^ /32 ||«||p 1 1) ^ 72 . (23) 
By choosing 1 < f3 2 < cr and 2 < r < 4, one can easily verify that -^z^h > 2 and (cr — 1)^^72 > 2. 
By (11) and (18), we obtain 

< ^IIM CT+1 ll 

= dMI^+i)- (24) 
Since 2(cr + 1) > 2 and (cr — 1)^3572 > 2 when cr > 1, by (12) we conclude 

/2<C a H|^||0i-«a|| r . (25) 
Now, when a = 1, by choosing 2 < r < 4 in (23) and applying inequality (24) with f3 2 = 1, we have 

Ili^OK- 1 !!^ = \\b{\v\^)\\^ 2 
< c a \\v\\i+Xi r 

By (12), we conclude that 

h<C a \\v\\%\\v x -v 2 \\ r . (26) 

Case 2. N > 3: 

By Holder's inequality, we have 



h < ||B(H <7+1 )||r'fl 1 |||«r 1 ||r'/J 1 ||«l-«2||r' 71 



II^M^JIUblip^JIt;!-^!!^, (27) 



where i + -f + i = 1. 

t>i pi 71 
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Now, by (8), (9) and (18), we have 

\\B(\v\° +1 )\\ r , ei < C\\B(\vr +1 )\\ w2 , m 
< C a \\\v\° +1 \\ m 

= C a,\\v\\^ +1)m , 

where pj- > ~ — when i — > 0, and r'6\ > m when ^ — < 0, and m > 1 to be decided later. 
Then we obtain 

h < CJv\\l^ 1) Jv\\l-\,^\v 1 -v 2 \\ r , J1 . (28) 
Now, by requiring (er + l)m = (a — l)r'0i = r'ji = r (choose j3\ = oo when a = 1), we have 

r - 1 

6*i = >l=^CT<r-l 

r — 1 — cr 

r - 1 

A = > 1 => cr < r 

(7 — 1 

71 = r — 1 > 1 
r 

771 = > 1 

(7 + 1 

N + 2 r - 1 

0~ < —^T7~ r — 1 or f > — ^ — • 

27V 2 
Since 2 < r < -^rj, we obtain cr < j^rg, i.e., 

l2<C a \\v\\ 2 ^\\v 1 -v 2 \\ r . 

The same can be shown for (21). □ 



Now, we are ready to prove Proposition 2. 
Proof. Recall that f(v) = B(\v\' 7+1 )\v\ cr ~ 1 v, by direct calculation, we have 

df(v) _ CT + 1 [B( | un|ur -i w + £ , (H .+i )|ur i ] 



dv 2 



when v ^ 0, where v* is the conjugate of the complex-valued function v. 
For v = 0, we have 



dv 
df(v) 



dv 

Now, by Lemma 3 and Lemma 4 we obtain 



0. 



("l - V 2 )\\ r > < C a \\v\\Hi\\vi -V 2 \\ r 



dv 

ll^^i-^Hr^CalHI^K-rall, 



dv 

When viv 2 7^ 0, by Mean- Value Theorem, we have 



\f(vi) - f(v 2 )\ < max{|^|, l^|}h - v 2 \. 
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That is, for some intermediate point vo between v\ and v 2 

||/(«l)-/(«2)||r < Ca||W0||^i||i;i-«2||r 

< C a (\\v 1 \\ H i + \\v 2 \\ Hl ) 2 ' J \\v 1 -v 2 \\ r . 
When Viv 2 = 0, without loss of generality, we assume that v 2 = 0, then 

||/(«l)-/(«a)||r' = \\m)-0\\ r > 

< C a (\\vi\\ H i + \\v2\\H^\\vi\\r, 

where the last inequality follows from direct application of Lemma 4. 
Therefore, we conclude that, for any Vi,v 2 G H 1 C L r , 

\\f(vi) ~ f(v2)\\r> < C a (II^H^ + \\v 2 \\m) 2a \\V1 - V 2 \\r ■ 

□ 

Next, we will give some elementary properties of the free evolution (linear Schrodinger equation) 
formally defined by the group of operators 

U{t) = exp(iiA), (31) 

where t 6 R. In the following, we will state some well-known results about the operator U(t) without 
proving them (see, e.g., [10], [21]). 

Lemma 5. For any r > 2, and for any t 7^ 0, U(t) is a bounded linear operator from L r to L r , and the 
map t — > U(t) is strongly continous. Moreover, for all t £B \ {0}, one has 

||^(t)»|| r <(47r|t|)^-^||t;|| tJ (32) 

for all v G L r . 

Corollary 6. Let I be an interval o/R 7 and let v £ C(J, L r ) . Then for allt S M the map r — > U {t—T)v{r) 
is continous from I\{t) into L r . 

3. Existence and Uniqueness of Local Solutions 

In this section we will prove a local existence and uniquenss theorem of solutions to system (4) by a 
fixed point technique. 
The integral equation 

v(t) = U(t - t )v Q + i [ U(t- T)f(v(r)) dr (33) 

J to 

may be considered as the integral version of the initial value problem for equation (4). 

Defining the subspace Y(I) C C(I,X) and Y b (I) C C b (I,X) by 

Y(I) ={«:i)6 C(I, X) and v(t) = U{t - s)v(s) for all s and t E 1} 
Y b (I) = Y(I)nC b (I,X). 

Here for special interest we choose the Banach space X = L r (IR iv ), for some r > 2, which is specified in 
the proof of Lemma 3 and Lemma 4, and X = L r (E^). 

If v G C b (I,X), we shall denote its norm by \v\i, and for v G C b (I,H l ), we denote its norm by \v\jj-l j. 
The ball of radius R in C b (I, X) will be denoted by B(I, R). 
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Let ti,t2 £ R and let v(t) be a family of complex- valued functions defined on R JV , depending on a 
parameter t 6 1. We formally define the operators 

[G(t 1 ,t 2 )v](t)=i [* U(t-r)f(v(T))dT, (34) 

where / is the nonlinear term defined in (13). The first lemma below gives a meaning to the expression 
defined by (34) and contains some of its properties. 

Lemma 7. For any interval I C M (possibly unbounded), the maps (tx, £2, v) — > G{t\, £2)1* are continuous 
from I x I x C(I,X) to Yb(K). Moreover, for any t±,t2 £ I, (ti < ^2), for any compact sub-interval J 
such that [£1,^2] C J C I , and for any t € [ti,t2]j for any v\,V2 G C{I,X) the G operator satisfies the 
estimates 

||G(tl,t2)l>l(*) - G(t U h)v2(t)\\r < k'\\vx - V 2 \\j\t2 - h\^^ +1 

where k' = fc(47r) ~ _ t , k = C a (\\vi\\ H i + I^Wh 1 ), which is derived in the proof of Proposition 2. 

Proof. For any v <E C(I,X) the function r —> f(v(r)) belongs to C(I,X) as consequence of Proposition 
2. Therefore, by Lemma 3, for any tgK \ {t} the function 

r^U(t-r)f(v(r)) (35) 

is continous from I to X. To check the intcgrability of the function (35) it will be enough to show the 
integrability of its norm. More generally one is interested in the integrability of 

\\U(t-T)\f(vi(r))-fMr))]\\r, (36) 
for any Vl ,V 2 € Ci^H 1 ) C C(I,X). 

This is a direct consequence of Proposition 2 and Lemma 3: For t£l, for every compact sub-interval 
J C I and r G J, we have 

\P(t - t)[/Mt)) - f(v 2 (T))]\\r> < (4rr|i - r|)£-**|h " «a|U ■ 
Finally, we come to the conclusion that 



||G(*i,t 2 )ui(i) - G(ti,ta)«2(t)|| < fc'K - «2||j|t 2 - ■ 



□ 



Now, in order to study the equation (33) one needs the operators 

[F(t )v](t) = [G(t ,t)v](t). (37) 
The existence and properties of F follow immediately from Lemma 5. 

For every v € C(I, H 1 ) C C(I, X), 

[A{t ,v )v}{t) = [F(t Q )v](t) + U(t - to)vo (38) 
is a continuous map from C(I, H 1 ) C C(I, X) into C(I, X). 

With these notations equation (33) may be rewritten as 

A(t ,v )v = v. (39) 

The next lemma gives some elementary properties of the solutions of equation (33). In particular, it 
expresses the consistency of the change of the initial time t n . 
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Lemma 8. Let I and J be two intervals of R, J C /, Zei to G J > let vq G if 1 be such that the function 
t —>U(t — to)vo belongs to Y(I), and let v G C(J,X) be a solution of the equation (39) 

(i) The function 

<f>(v):s->U(--s)v(s) = [<i>(v)](s) (40) 
belongs to C(J,Y(I)) and satisfies for all s,s' G J the equality 

lcf>(v)}(s)-[cf>(v)}(s') = G(s',s)v. (41) 

Furthermore, if for some s G J, [(f>(v)](s) G Yf,(I), then (j)(v) G C(J, Yj,(/)). If in addition J is bounded, 
then 4>{v) G C b (J,Y b (I)). 

(ii) For any s G J, u satisfies the equation 

A(s,v(s))v = v. (42) 

Proof. Apply the operator U(t — s) to equation [A(to, Vo)v](s) = v(s) and use the fact that U(t — 
s)[G(ti,t2)v](s) = [G(ti,t2)v](t) (for the proof of this identity, we refer to Ginibre and Velo [10]) yields 

[#»)](*) = U(- - t )vo + G(t , s)v. (43) 

^From which (41) follows immediately. The continuity properties of the left-hand side of (43) are then a 
consequence of the assumptions made on vq and of Lemma 7. Finally, putting s' = t in (41) and taking 
the values of both members at t one obtains equation (42) at time t. □ 

We are now ready to discuss the problem of the existence and uniqueness of solutions of equation (39). 

Theorem 9. For any p > 0, there exists a Tq(p) > 0, depending only on p, such that for any to G K 
and for any Vq G H , for which \\vq\\hi < p, equation (39) has a unique solution on C(I,X), where 
I=[h- T (p),t + T (p)] and X = U. 

Proof. Let p be a fixed postive number, let t and T G R, T > 0, and let I = [to — T, t + T]. Then for 
every Vi,V2 G H 1 and — «2||hi < 2p, Lemma 5 and (37) yield the inequality 

\F(t )vi - F(t )v 2 \i < 2k'\t - to|"-^ + >i -v 2 \i. (44) 
In particular, if we take T = To(p) with 7o(p) defined by 

4fc'|T (p)-io|" +1 = l (45) 

in equality (44) it gives 

\F{t )vi - F{t Q )v 2 \i < -\v! - v 2 \i . (46) 
Let now vq G X be such that U(- — to)vo € B(I,p). Definition (38) and estimate (46) imply 

\A(to,vo)v\i<2p, (47) 

and 

\A(t ,v )vi ~ A(t ,v a )v 2 \i < hv! - v 2 \i , (48) 

for all v, v\,V2 G B(I, 2p), from which it follows that A(to,vo) is a contraction from the ball B(I, 2p) into 
itself. The result is now a consequence of the contraction mapping theorem. □ 
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4. Global Existence of Solutions 



In this section we will study global existence of solutions to system (4) under the condition of a > 1. 
We will show below that we have global solutions when 1 < a < -4. Comparing this to the results of 
the classical NLS (a < j^), we "gain" global regularity for larger range of values of a. As we stated in 
the introduction, system (4) will recover the classical NLS as the parameter a — > + , in a subscquencial 
paper, we will study numerically the blow up profile of the classical NLS by focusing on system (4) with 

< a < 4f when N < 3. To be more specific, the profile of blow-up in the critical case a = jj in the 
classical NLS has not been known completely, in a subsequcntial work, we will compute SH system (4) 
and try to find out the blow up profile by forcing the parameter a to approach zero. 

Theorem 10. Let vq 6 H (R ). If 1 < a < -4, then there exists a unique solution 

v G C((— oo, oo); 7? 1 (IR Ar )) of the initial-value problem (4), in the sense of the equivalent integral equation. 

Furthermore, as long as v(x,t) remains in H 1 (M. N ) ; the energy 



Af(v) 



\v(x, t)\ 2 dx 



(49) 



and Hamiltonian 

H(v) 

remain constant in time. 



2 U{x,t)\v{x,t)\ a+1 



(7+1 



dx 



(50) 



In the local existence theorem in section 3, we have shown that the length Tq, of the interval of 
existence [io,io + To], can be taken to depend only on ||uo||ifi. It follows that if v(x,t) is a maximally 
defined solution on [to,Tmax)j then either 

Tmax = +oo 



or 



lim 1 1 f (t) 1 1 jyi = +oo. 



max 



The heart of the global existence proof lies in the use of the invariants (49) and (50), which enable us to 
obtain an a priori bound of the following type: 



\\v(x,t)\\ H , <C{N,H). 



(51) 



Proof. We proceed as follows: 
^From (50) we have 



Observe that 



\\vv{x,t)\\ 2 <n 



(7+1 



u\v\ a+i dx 



(52) 



where 



r +i dx 



< IMUK +1 II P ' 



h\ p \m;% +1) 



- + ^- = 1, 1 < p, p' < OO. 

p pi 



(53) 
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Case 2. N < 2: 

By (8) and (18), we have 



\\u\\ p < C\\u\\ W 2, P 

= c\\B(\vr +i ) 



< C a \\\v\' 



tr+1 
CT+1 



for any p > 1 . 

Plug the above inequality into (53) and requiring p = p 1 = 2, we obtain 



By (7), taking q = 2(a + 1), we obtain 



< C a \\v 



|2(<T+1) 

2(cr+l) 



l«l| 3 (a+i) < C||«irwtTr w ||v w ||wtn 



with 



o< 2( ^ +1) - 2 iv<i ; 



2-2(ct+1) 



which is always satisfied when N < 2. 
Then (52) yields 



||Vi;(t)|| 2 <H a + C a \\v \\ 2 ^-^ N -V\\Vv\\ 



aN-2 



Case 2. N > 3: 

By (9) and (18), we obtain 



where 



\\u\\ P < C\\u\\ W 2, m 

= C\\B(\v\^)\\ w .,. 

< C a \\\v\° +1 \\ m 

= Coc\\v\\\v +1 ) m , 



112 N 

~ = "77 > ^ m < TT • 

p m Jy 2 



Pluging into (53) and requiring m = p' , i.e., m = j^rg, we get 



2((T+1) 



By (7), taking q — m(a + 1), wc obtain 



i(°-+l)-2 



+ 1) ' 



-AT 



with 



Then (52) yields 



„ m(cr + l)-2 Ar „ 4 
< — '-^N <!=> a < 



2m((j + 1) 



N -2 



||Vv(t)|| 2 < Ho + C Q || Vo || 2 (ff+1) - (CTAr - 2) ||V^| 



aN-2 
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For (54) and (57), ||Vu|| is bounded when aN - 2 < 2, i.e., a < ±. Therefore, the H 1 norm of 
the solution v is bounded uniformly independent of time t, so we can conclude that we have global 
solution for any 1 < a < . □ 

In conclusion, we have shown that the Schrodinger-Newton system (3) and the Schrodingcr-Hclmholtz 
system (4) admit short time unique solution when 1 < a < N 4 _ 2 (by definition, N A _ 2 = oo when N < 2) 
and global existence of unique solution when 1 < a < -4. Comparing to the result of classical NLS (1) 
(<t < jf^2 f° r l° ca l existence and a < for global existence), one expects this "better" result since the 
nonlinear terms in system (3) and system (4) are milder than that of the classical nonlinear Schrodinger 
equation (1). 
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